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ABSTRACT 


■ — The  set  partitioning  model  is  used  as  the  basis  for  an  interactive 
approach  for  solving  a  broad  class  of  routing  problems.  A  pricing  mecha¬ 
nism  is  developed  which  can  be  used  with  a  variety  of  methods  in  generat¬ 
ing  improving  solutions.  A  version  of  the  approach  has  been  implemented 
via  a  colorgraphics  display. 
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1.  INTRODUCTION 


We  will  consider  here  a  set  partitioning  based  approach  for 
soLving  a  broad  class  of  routing  problems.  The  approach  is  designed 
to  take  advantage  of  a  high  level  of  human  interaction;  the 
current  implementation  is  interactive  via  a  colorgraphics  display. 
However,  many  of  the  concepts  discussed  here  could  be  easily  implemented 
in  an  automatic  system. 

The  routing  problem  which  motivated  much  of  this  work  is  what 
is  called  the  static  or  subscriber  dial-a-ride  problem.  This  problem 
will  be  discussed  in  detail  in  later  sections.  It  is  one  of  the  more 
complex  members  of  the  class  of  routing ' problems  which  are  amenable  to 
the  approach  presented  here.  'This  class  also  includes  many  practical, 
delivery  problems.  In  order  to  introduce  the  underlying  methodology 
which  provides  the  basis  for  the  approach,  consider  a  very 
simple  delivery  example.  Assume  that  a  depot  is  located  at  the 
square  box  labeLed  D  in  Figure  1.  From  this  depot  a  single  delivery 
is  to  be  made  to  each  of  the  points  represented  by  numbered  circles. 

The  numbers  on  arcs  connecting  the  circles  represents  the  travel  dis¬ 
tance  between  delivery  points.  Assume  also  that  each  vehicle  (e.g. 

truck)  can  deliver  to  a  maximum  of  two  points  on  a  single  trip.  The 
objective  is  to  determine  which  vehicle  should  deliver  to  each  point 
and  the  routing  for  the  vehicles  which  minimizes  the  total  distance 
t  ravi  1  led  . 

F.ach  column  in  the  matrix  o[  Table  1  represents  one  possible  vehicle 
route.  For  example,  column  one  represents  a  vehicle  travelling  from  the 
depot  to  delivery  point  (1)  and  returning.  The  c ,  row  indicates  the 
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distance  traveled  for  each  trip.  For  example,  column  six  represents 
a  vehicle  proceed  in;;  from  the  depot  to  delivery  point  (1),  on  to  delivery 
point  (2)  and  from  there,  back  to  the  depot.  The  value  of  c^  is  14, 
the  total  distance  traveled  for  this  trip.  By  enumerating  each  of  the 
possibilities,  as  has  been  done  for  this  matrix,  the  problem  becomes 
one  of  selecting  a  set  of  columns  such  that  every  row  is  represented 
in  exactly  one  column  and  the  sum  of  the  costs  of  the  columns  selected 
is  the  smallest  possible.  This  integer  program  is  called  a  "set  parti¬ 
tioning  model". 

The  set  partitioning  model  was  originally  proposed  for  the  delivery 
problem  by  Balinski  and  Quandt  [2].  The  model  is  very  powerful  in  the 
sense  that  many  realistic  route  constraints  and  cost  functions  can  be 
handled  easily  in  the  column  enumeration  process.  The  obvious  short¬ 
coming  of  the  model  is  that  there  are  typically  a  very  large  number  of 
columns  to  be  enumerated  and  the  resulting  integer  program  is  very 
large.  The  approach  presented  here  is  heuristic  in  the  sense  that  we 
generate  only  a  subset  of  the  possible  columns  or  routes  and  in  general 

we  do  not  solve  the  set  partitioning  model  to  optimality. 

The  set  partitioning  model  has  two  very  desirable  features  for 

interactive  optimization.  The  first  is  that  any  route  generated  can  be 
included  as  a  column  in  the  model.  This  allows  the  human  interactor 
to  ul l It/e  his/her  Intuition  and  spatial  preceptlon  as  well  as  a  wide 
spectrum  of  mathematical  techniques  to  generate  new  routes.  The  second 
feature  is  that,  unlike  more  general  integer  programs,  a  feasible 
solution  to  the  set  partitioning  model  provides  the  basis  for  pricing 


information  which  can  be  used  to  generate  new  candidate  columns. 

We  will  restrict  the  class  of  routing  problems  considered  here 
to  be  those  for  which  any  subroute  of  a  feasible  route  is  also  a 
feasible  route  with  cost  less  than  or  equal  to  the  cost  of  the  original 
route.  By  imposing  this  restriction,  as  long  as  there  is  at  least  one 
1  in  every  row  of  the  partial  set  partitioning  model  that  we  have 
enumerated,  we  can  easily  generate  feasible  partitions.  The  only 
other  restriction  that  we  put  on  the  class  of  routing  problems  is 
that  we  be  able  to  pose  them  in  a  natural  way  as  set  partitioning 
problems.  However,  it  should  be  noted  that  if  there  is  not  a  nice 
spatial  representation  of  the  routing  problem,  the  human  interactor 
is  much  more  restricted  in  his/her  contribution. 
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2.  SET  PARTITIONING  AND  ROW  PRICES 


The  set  partitioning  problem  can  be  stated  as 


minimize 


=  y  c.x. 

A  J  J 


subject  to  y  a  .x.  =  1  for  i  =  1,2,... ,m 

j=l  J  3 

x.  =  Oor  1  for  j  =  l,2,...,n. 


For  the  delivery  example  in  Table  1,  the  set  partitioning 


model  has  the  second  row  as  the  value  of  the  c.  and  rows  three 

3 

through  seven  as  the  values  of  the  a.  ^ .  The  variables  which  are 

set  to  one  in  a  solution  to  the  set  partitioning  problem  will  be 

k  k 

called  a  "partition."  We  will  denote  a  partition  as  J  =  {j|x.  =  1}. 
Balas  and  Padberg  [1]  provide  a  recent  survey  of  results  related  to 
set  partitioning  problems. 

A  fundamental  idea  underlying  much  of  the  work  presented  here 
is  the  concept  of  "row  prices." 

Definition:  P"^  =  (p^p^, . . .  ,p^)  is  a  set  of  feasible  row  prices 
corresponding  to  the  partition  if 


It  will  be  useful  to  interpret  the  price  p^  as  an  estimate  of  the 
cost  to  satisfy  constraint  i  using  solution  X^.  For  the  delivery  problem, 
p^  is  then  an  estimate  of  the  cost  of  satisfying  the  requirement  of 
delivery  point  i  using  the  route  corresponding  to  partition  J’*'. 

Theorem  1:  Given  a  set  of  feasible  row  prices  (p^ ,p^ , • • •  , p^)  corresponding 
to  partition  J  with  value  z\  any  other  partition  has  value 


■  zl  -  ;j2 

jeJ 


l 

i=l 


(p . a . ,-c.) 
J-  il  J 


Proof : 


ll 

j£j 


<  1 

i  -1 


P. a. ,-c.) 

1  ij  3 


m 

1 

i=l 


u 

jeJ 


iJ 


2  r 

Since  J  is  a  partition,  l  a  -  1  for  each  i  =  l,2,...,m.  Also, 

jeJ  ^ 

since  (p^.p^, . . . ,p^)  are  feasible  row  prices  corresponding  to  X1  we 

r  1  1 

have  J  p.  =  Z  .  Hence  the  result  follows. 
i=l  1 

Corrol lary  1:  For  any  set  of  feasible  row  prices  P1  corresponding 

1^1  1 
to  a  partition  J  if  \  (p . a . . -c  )  ^  0  for  j  =  1,2, ...,n  then  X  is 

i=l  1  1J  J 

optimum . 

It  can  also  be  shown,  using  linear  programming  duality,  that  a 
set  of  feasible  row  prices  P^1  satisfying  Corrollary  1  exists  if  and 
only  if  X1  is  an  optimum  solution  with  the  constraints  =  0  or  1 
replaced  by  x^  >_  0  for  j  =  l,2,...,n. 


•  ■  ,,£<  y'&.+r.  , 


■;  .  V 
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The  quantity  J  p ~ a . .  -  c.  will  be  interpreted  as  the  "potential" 
i=1  1 

savings  over  the  value  of  Z  which  can  result  from  constructing  a  parti¬ 
tion  that  includes  column  j.  Note  from  Theorem  1  that  the  potential 


savings  can  actually  be  achieved  only  if  a  partition  can  be  constructed 


from  columns  with  nonnegative  potential  savings. 
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POTENTIAL  SAVINGS  HEURISTIC 


Given  a  partition  1  a  corresponding  set  of  feasible  row  prices 
p\  an  attractive  heuristic  for  attempting  to  generate  a  better  parti¬ 
tion  is  the  following: 

2  2 

Step  (0):  Let  J  =  0  (J  will  be  the  indices  of  columns  in  the  new 
partition)  and  N  =  {l,2,...,n},  (N  will  be  the  indices  of 


columns  which  are  candidate  for  inclusion  in  J  ) 

m 

Step  (1):  Calculate  the  potential  savings  J  p^a..  -  c.for 


i=l 


ij 


j  =  1 , 2 , . . . ,  n . 


Step  (2) :  Determ 


mi  m 

r  I  P  I 

ine  keN  such  that  >  p.a.,  -  c  >  i  pa 

i  u  — 


i=x  '  "  “  i=l  1  J 

For  all  jeN  (i.e.,  pick  the  column  in  N  with  the  largest 


potential  savings) 

Step  (i):  11  1  2ai]  =  chen  set  =  0.  (Note  from  the  assumption 

of  section  1  that  anv  subroute  of  a  feasible  route  is  also 


Step  (4) 


a  feasible  route  the  new  column  is  legitimate) 

2  2 

Let  J  =  .1  IJ  f  k }  (i.e.,  put  column  k  in  the  new  partition) 


Step  (5):  Delete  from  N  all  j  for  which  a., 

lk 


1  and  a . .  =1  for  some 
1J 


i  =  1,2, 


>  m . 


Step  (b):  If  N  =  0  stop.  Otherwise  go  to  step  (2). 

Note  that  under  the  assumption  that  any  subset  of  a  route  i6 

also  a  feasible  route  (discussed  in  section  1)  this  procedure  will 

2 

always  terminate  with  J  as  a  partition  although  not  necessarily  a 
better  partition  than  J  .  In  addition,  computation  to  date  indicates 
that  an  optimum  or  near  optimum  solution  to  the  set  partitioning 
problem  is  determined  verv  quickly  by  repeated  application  of  the 
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potential  savings  heuristic.  The  heuristic  is  repeated 

j, 

until  either  optimality  is  proven  (i.e.,  for  some  X  all 
m  k 

I  Pj  a,  .  -  c.  <_  0)  or  until  some  specified  number  of  partitions 
i  =  l  J  J 

has  been  generated. 

To  illustrate  the  potential  savings  hueristic,  consider  again 

the  delivery  example  depicted  in  Figure  1  and  Table  1.  Suppose 
that  we  select  =  {l, 2,..., 5}  as  an  initial  partition.  A  set  of 

feasible  row  prices  P  *  is  given  in  Table  2.  (The  question  of  how 
to  generate  "good"  feasible  row  prices  will  be  addressed  in  the 

<»  p 

next  section.)  The  corresponding  potential  savings  J  p‘.a  .  -  c. 

•■*  1  i  i  i  -i 


i=i 


i  ij 


are  also  given  in  Table  2.  Applying  the  potential  savings  heuristic 


and  breaking  ties  by  selecting  the  column  with  the  lowest  index 

2 

yields  the  new  partition  J  =  {6,13,5}.  The  new  partition  has  a  cost 

of  =  22  as  compared  to  a  cost  Z*  =  28  for  the  initial  solution. 

2  ?  2 

Using  the  row  prices  P  and  potential  savings  l  p^a^,  -  c. 

L=1  J  J 

of  Table  2  and  reapplying  the  potential  savings  heuristic  yields  the 

3  3 

partition  .1  =  {8,10,5}  which  has  a  cost  of  Z  =  20.  Again,  from 

3 

Table  2  we  find  that  using  the  row  prices  P  given  potential  savings 
m  i 

£  p  -i .  .  -  c.  <0  for  j  =  l,2,...,n.  Hence,  from  Corroilary  1  the 
1=1  .  .  3 

partition  .1  is  optimum. 
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4.  ROW  PRICING 


For  a  given  partition  X  a  set  of  feasible  row  prices  is  obtained 

by  allocating  the  column  cost  c  for  each  jeJ  among  the  rows  having 

j 

a^j  =  1.  For  the  delivery  example,  this  corresponds  to  allocating 

the  trip  cost  among  the  delivery  points  of  the  trip.  When  a  column 

k  k 

j e J  dontains  only  one  a..  =  1,  the  row  price  is  p .  =  c..  However, 

ij  1  J 

when  a  column j  contains  more  than  one  a..  =1,  there  are  an  infinite 
number  of  possible  sets  of  prices.  As  an  example  consider  the  parti¬ 
tion  J3  =  {8,10,5}  for  the  problem  in  Table  2.  Since  column  5  has 

only  a.  „  =  1  the  value  p^  =  2  is  unique.  Column  8  has  both  a,  0  =  1 

J  >  J  J  lyO 

and  a.  Q  =  1,  hence  the  cost  c.  =  8  could  be  allocated  between  rows 

4 ,  o  o 

1  and  4  in  an  infinite  number  of  ways.  Similiarly,  column  10  has 
both  a2  =  1  and  a^  =  1,  hence  c^q  =  10  could  be  allocated 

between  rows  2  and  3  in  an  infinite  number  of  ways.  If  we  allocate 

3  3  3  3 

Cg  as  p^  =  4  and  p^  =  4  and  allocate  c^0  as  p^  =  5  and  p^  =  5,  the 

m  3  3 

resulting  J  p.a.  -  c.  do  not  indicate  that  J  is  an  optimum  parti¬ 
al  3  3  3 

tion.  Hence  the  set  of  prices  P  given  in  Table  2  are  clearly  better 

3 

since  they  do  indicate  that  J  is  on  optimum  partition. 

Ideally,  we  would  like  a  set  of  prices  which  would  drive  the 
potential  savings  heuristic  toward  an  improving  solution  and  would 
indicate  optimality  when  no  improving  solution  is  possible  (i.e.. 

We  would  like  the  prices  to  be  analogous  to  dual  variables  in  linear 


programming).  Unfortunately,  it  is  easy  to  construct  cases  for  which 


no  such  prices  exist  (i.e.,  any  problem  for  which  the  integer  and 
continuous  solution  differ).  For  the  delivery  problem  and  the  more 
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5.  COLUMN  CKNfc! RATION 

1  For  large  scheduling  and  routing  problems  it  is  generally  not 

practical  to  generate  all  columns  of  the  corresponding  set  partitioning 
model.  The  remainder  of  this  paper  will  be  concerned  with  using  in¬ 
formation  gleaned  from  one  soLution  via  Theorem  1  to  generate  a  new 
and  hopefully  better  solution.  This  is  accomplished  bv  either  generating 
new  columns,  adding  them  to  the  current  set  partitioning  model,  and 

*  then  resolving  the  model  or  by  using  the  information  from  Theorem  1 
directly  to  generate  a  new  solution.  In  the  latter  case  it  is  not 
necessary  to  retain  the  columns  of  the  set  partitioning  model.  How- 

*  ever,  if  the  columns  are  retained,  it  is  possible  to  further  improve 
the  solution  by  periodically  solving  the  set  partitioning  model. 

We  should  note  that  for  the  class  of  scheduling  and  routing  prob- 
^  lems  being  considered  here,  it  is  very  easy  to  generate  an  initial 

solution.  For  our  examples  we  use  the  identity  solution  (e.g.,  in 
the  delivery  problem  this  is  the  solution  which  has  each  vehicle  making 

*  a  single  delivery)  as  our  initial  solution.  However,  any  feasible 
solution  could  be  used  as  the  initial  solution. 

Clearly  there  is  a  broad  spectrum  of  possible  approaches  that 

^  one  might  use  to  generate  new  columns  and/or  new  solutions  to  the 

set  partitioning  model.  In  fact,  variations  of  many  of  the  heuristics 

which  have  been  applied  to  delivery  problems  can  be  used  very  effectively 

^  in  conjunction  with  Theorem  1.  In  the  next  section  we  will  discuss 

the  use  of  the  Clarke  and  Wright  [3]  savings  procedure  in  conjunction 
with  the  delivery  problem.  In  later  sections  we  will  discuss  more 

|  complex  clustering  and  chaining  heuristics  as  we  have  applied  them 

to  the  diai-a-ride  problem. 


Cl.ARKh  AND  WHICH  I  l  '.vi  lli  PRICINC 


tv . 


Ilu  "s.ivinKs"  heuristic  ui  Clarke  . 1 1 1 J  Wright  |  J  )  is  tin1  most 
widely  known  ot  tin1  heuristics  developed  to  date  tor  delivery  problems. 
Ihe  algorithm  proceeds  bv  calculating  a  savings  for  each  pair  of  deli¬ 
very  so  inis  i  and  |  di  lined  as 


(d,,.fd.  -HI  ) 

()i  ij  j(l 


d  . 
i  1 


wit  it  It  is  t  be  savings  in  mileage  ot  supplying  delivery  points  i  and  j 

on  tin1  same  route  as  opposed  to  supplying  them  individual lv  directly 

t  io.ii  l  be  depot  td|(.  is  the  distance  I  rotn  the  depot  to  delivery  point  i). 

Routes  are  then  eonsirueted  either  one-a t -a -t ime  or  in  parallel  bv 

eons i.ier i ng  pairs  of  points  in  order  ot  deereasing  savings  and  including 

them  in  the  same  route  it  sueh  a  route  is  leasihle. 

Suppose  that  we  eoiisider  ouee  more  the  delivery  example  in  Figure  1 

and  the  covering  solution  information  in  Table  ' .  Note  that  the  Clarke 

and  Wright  savings  values  are  exactly  the  values  of  /  pa  -  c  in 

,  i  =  l  1  U  1 

fable  d  since  d  .  -  p  .  /  d  tor  i  -  l,d,...,m.  Applying  t  lie  C'-W  savings 

algoiithm  yields  the  same  routing,  eonf  i  gura  t  i  on  as  d~  =  At 

this  point,  the  C-W  algoi i t hm  would  terminate.  However,  suppose  that 

we  set  d().  =  Pj/2  for  i  =  l,d,...,m.  The  new  savings  are  then  exactly 
, .  > 

the  /  pi,.  -  e,  in  Table  d.  Applying  C-W  algoiithm  yields  the 
1=1  1  1 !  1 

same  routing  configuration  as  I  -  {H,l(),si  which  as  noted  previously  is 
the  optimum  solution  to  this  delivery  example.  Hence,  for  this  example 
at  least  the  C-W  algorithm  without  pricing  d  id  not  vield  an  optimum  solution 
while  the  same  algorithm  when  combined  with  pricing  did  vield  the  optimum 


solut ion. 


Note  that  when  routes  are  allowed  to  contain  more  than  two  trips, 
onlv  a  small  subset  of  the  set  partitioning  columns  would  be  generated 
using  this  procedure.  Also,  note  that  the  prices  and  savings  can  be 
calculated  without  ever  generating  the  set  partitioning  matrix. 

When  the  number  of  points  allowed  in  a  route  exceeds  two,  some 

interesting  questions  arise  as  to  exactly  how  the  algorithm  should  be 

implemented.  As  an  illustration,  suppose  that  in  the  example  we  allow  a 

vehicle  to  deliver  to  at  most  three  points  rather  than  two.  Now  suppose 

3 

that  we  start  with  the  solution  J  in  Table  2.  We  see  that  the  two-at-a- 
t ime  potential  savings  ate  all  non-positive.  However,  if  we  ignore  this 
fact  and  proceed  with  the  algorithm,  we  would  put  points  1  and  4  in  the 
same  route  since  their  potential  savings  of  zero  is  maximum.  If  we  then 
consider  adding  a  third  point  to  this  route,  we  find  that  adding  point  5 
has  a  potential  savings  of  2  which  is  maximum.  Finally,  we  combine  points 
2  and  1  into  a  single  toute  since  this  has  a  potential  savings  of  zero. 

The  resulting  routes  (4,15)  and  (2,3)  has  a  length  of  10.  Therefore, 
it  appears  that  we  can  get  better  information  by  recalculating  the  poten¬ 
tial  savings  after  each  augmentation  of  a  route.  This  recalculat  on  is 
not  done  in  most  implementations  of  the  C-W  algorithm. 

Clearly,  when  constructing  a  route  containing  more  than  two  points 
one  must  decide  where  in  Lire  route  to  put  each  additional  point.  The  poten 
tial  savings  can  he  determined  exactly  only  by  solving  a  travelling  sales¬ 
men  problem  over  each  new  point  which  is  a  candidate  to  be  added  to  the 
route.  This;  is  computationally  expensive  if  a  route  can  contain  a  large 
number  of  points.  Tn  most  implementations  of  the  C-W  algoritiim,  new  points 
are  simply  added  on  to  the  end  of  the  route  being  constructed.  When  the 


f 


» 


» 


) 


* 


algor  itlun  is  imp lomentej  interne Lively  using  computer  graphics,  it  appear 
that  the  human  can  perforin  an  important  role  both  in  selecting  candidate 
points  and  in  inserting  them  logically  into  routes. 
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7.  LOCATION  -  AL LOCATION  WITH  PRICING 

Another  procedure  which  Krolak  and  Nelson  [5]  have  found  effective 
in  approaching  the  delivery  problem  utilizes  the  location  -  allocation 
model  ol  Cooper  [4],  This  basic  concept  can  also  be  used  in  conjunction 
with  Theorem  1  to  generate  intuitively  appealing  columns  to  add  to  the 
set  partitioning  model. 

To  illustrate  the  procedure  consider  tire  delivery  example  illustrated 
in  figure  J  .  Again  the  circles  represent  delivery  points  and  the  square 
represents  the  depot.  The  basic,  idea  is  to  use  a  surrogate  distance 
rather  than  the  actual  distance  in  determining  the  points  which  are  assigned 
to  each  vehicle.  Hie  surrogate  distance  is  obtained  by  assuming  that  the 
vehicle  travels  from  the  depot  to  a  specified  cluster  point,  represented 
in  Figure  2  bv  Lite  dashed  circles.  It  then  makes  the  deliveries,  returning 
after  each  delivery  to  the  cluster  point.  After  all  deliveries  have  been 
made,  tin  vehicle  returns  to  the  depot.  Under  this  surrogate  distance, 
the  piohlem  becomes  one  ol  locating  t  lie  cluster  points,  one  for  each 
vehii  i e ,  and  then  assigning  the  delivery  points  to  each  vehicle. 

It  (a(|,l>())  represents  the  coordinate',  of  the  depot,  represents 

the  coordinate-,  oi  delivery  point  i,  (x  ,v.)  represents  the  coordinates 
ol  i' luster  point  i,  and  assuming  Euclidean  distance,  the  problem  can  be 
tnode  1  ed  as  follows: 


™  n  ,  ,  i 

min  /  /  d  |  ( x  .  -a  )  +  (y  . -b  )  ^  J  * /.  .  . 
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where  n  is  Che  number  of  vehicles  and  K  is  the  vehicle  capacity.  When 
z.  .  =  1,  delivery  point  i  is  assigned  to  vehicle  j  and  when  z  =0, 

ij 

delivery  point  i  is  not  assigned  to  vehicle  j. 

While  there  is  no  method  for  efficiently  solving  this  problem 

optimally,  the  following  is  an  attractive  heuristic.  Pick  a  set  of 

locations  for  the  cluster  points.  With  these  coordinates  fixed,  solve 

the  resulting  assignment  problem.  With  these  values  of  z  .  fixed, 

ij 

solve  the  resulting  location  problem.  Continue  alternating  between 
the  assignment  and  location  problems  for  some  specified  number  of 
iterations  or  until  there  is  no  further  improvement  in  the  objective. 
Once  a  cluster  has  been  determined,  the  vehicle  is  then  routed  among 
the  points  of  the  cluster. 

A  slight  modification  of  this  model,  together  with  Theorem  1, 

allows  us  to  generate  attractive  new  columns  for  the  set  partitioning 

problem.  Suppose  that  we  have  a  solution  to  the  set  partitioning 

problem  and  a  set  of  row  prices  p  ,p  ,...,p  .  Now  consider  the  model 

1  z  m 


t 


21 


E  r  I 

f 

i;  i 

:  t 

I 

!  f 


m 


I 

i-l 


l  P  • z  •  ■ 

j-1  1  1J 


s.  t . 


1=1  1J  - 


<  K  j  =  1,2, 


,n 


n 

y  Z. .  <  1  i  =  1,2, 

A lj  - 


,m 


z .  .  =  0  or  1 . 
ij 


I  ’ 

i 

I 


/ 

i 

t 


l 


*  » 

f 


I 


Any  cluster  generated  by  this  model  will  have  a  positive  potential 
savings,  with  respect  to  the  surrogate  distance.  Hence,  it  corresponds 
to  an  attractive  column  to  add  to  the  set  partitioning  problem  (consider¬ 
ing  surrogate  distances). 

Since  the  second  constraint  has  been  changed  to  an  inquality,  not 
all  delivery  points  will  be  assigned  to  cluster  points.  This  simply 
means  that  the  current  row  price  for  the  delivery  point  is  more  attrac¬ 
tive  than  the  cost  of  serving  the  delivery  point  in  alternative  clusters 
considered  by  the  model.  The  model  can  be  solved  using  the  same  heuristic 
discussed  for  the  earlier  location  -  allocation  model. 
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8.  DIAI-A-RIDE  PROBLEM 
i 

The  dial-a-ride  problem  is  a  much  more  complex  routing  problem 
than  the  delivery  problem.  In  the  dial-a-ride  problem  we  are  given 
an  origin-destination  trip  matrix  and  an  underlying  network  on  which 
the  trips  are  to  be  made.  There  is  a  single  item  (people,  goods, 
etc.)  (demand  for  service)  at  each  origin  that  needs  to  be  transported 
to  its  specitied  destination.  The  items  are  transported  from  origins 
to  destinations  on  vehicles  each  having  capacity  K.  We  wish  to  satisfy 
the  trip  requirements  while  travelling  the  minimum  distance. 

,  This  is  a  "static"  version  of  the  dial-a-ride  problem  since  time 

if,  not  considered .  There  are  a  number  of  more  comolex  versions  of 
ti;is  problem,  but  this  version  is  sufficient  to  demonstrate  the  basic  ideas 
j  of  our  approach. 

The  set  partitioning  model  for  the  dial-a-ride  problem  is  analogus 
to  that  of  the  delivery  problem,  but  here  rows  represent  trips  rather 
I  than  simply  delivery  points.  The  vehicle  capacity  constraints  are  handled 

by  generating  only  routes  which  satisfy  them. 

> 

» 
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9.  DECOMPOSITION 


In  delivery  problems,  representative  of  the  "one-ended"  class  of 
routing  and  scheduling  problems,  we  need  only  be  concerned  with  a  single 
stop  for  a  vehicle  to  satisfy  a  particular  demand  for  service.  In  con¬ 
trast,  the  "two-ended"  class,  which  includes  dial-a-ride  problems,  requires 
two  stops  in  a  specific  order  (sequence)  to  satisfy  a  specific  demand  for 
service.  It  is  this  requirement  for  sequencing  of  "pickup"  and  "dropoff" 
point  pairs  which  adds  greatly  to  the  difficulty  of  handling  the  two-ended 
class  of  vehicle  routing  problems. 

When  one  considers  examples  of  two-ended  problems,  it  becomes  imme¬ 
diately  apparent  that  the  sequencing  requirements  greatly  inhibit  the 
complex  pattern  processing  abilities  of  the  human.  Figure  3  illustrates 
an  example  of  a  25  trip  dial-a-ride  problem.  Rather  than  displaying 
order  and  structure,  the  problem  resembles  so  much  spaghetti .  It  is 
clear  that  lor  such  problems  the  human  intorac Lor  needs  more  help  in 
generating  good  columns  lor  the  set  partitioning  model. 

Unfortunately,  it  is  also  more  difficult  to  apply  straight  forward 
methods  such  as  the  savings  approach  discussed  earlier  for  the  delivery 
example.  In  generating  a  dial-a-ride  route  one  must  be  concerned  with 
where  both  the  origin  and  the  destination  occur  in  the  sequence  in  order 
to  calculate  the  potential  savings.  In  addition,  the  capacity  constraint 
may  negate  what  otherwise  appears  to  be  good  positions  for  the  origin 
and  destination  in  the  sequence. 

Because  of  this  complexitv,  it  is  helpiul  to  "decompose" 
the  problem  into  two  levels  which  we  call  "clustering"  and  "chaining". 


In  essence,  we  consider  a  route  to  be  made  up  of  two  components. 


Clusters 
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correspond  to  trips  which  can  all  be  on  a  vehicle  at  one  time,  while 
chains  correspond  to  movement  from  the  end  of  one  cluster  to  the  beginning 
of  the  next.  Figure  4  provides  an  example  of  five  good  clusters,  while 
Figure  5  illustrates  one  way  in  which  these  five  clusters  might  be  linked 
into  two  chains. 

The  partitioning  model  and  pricing  concepts  can  be  effectively  exploited 
to  aid  in  generating  improving  clusters  and  chains  in  a  column  generation 
approach  to  solving  two-ended  vehicle  routing  problems.  One  partitioning 
model  can  be  utilized  in  generating  improving  clusters  while  another  parti¬ 
tioning  model  helps  identify  better  chains. 

In  the  partitioning  model  for  clustering,  the  columns  represent  clusters 
and  the  rows  represent  the  individual  trips  (demand  for  service).  We  shall 
demonstrate,  in  the  next  section,  just  how  the  pricing  information  from  the 
partitioning  problem  can  be  used  to  generate  additional  clusters. 

In  the  partitioning  model  for  chaining,  the  columns  represent  chains 
and  the  rows  represent  the  individual  trips.  The  function  of  the  chaining 
partitioning  model  is  to  combine  the  clusters  into  good  vehicle  routes. 

The  overall  solution  procedure  consists  of  linking  the  clustering 
models  and  chaining  models  together  in  an  interactive  manner.  The 
clustering  models,  partitioning  matrix  and  pricing  information  are  used 
to  generate  good  clusters.  These  clusters  are  then  passed  to  the  chain¬ 
ing  phase  where  they  are  linked  together  via  the  chaining  models,  par¬ 
titioning  matrix  and  associated  pricing  information.  After  chaining, 
it  is  possible  to  return  to  the  clustering  phase  to  identify  additional 
clusters . 

The  next  two  sections  discuss  the  specifics  of  clustering  and  chaining 
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10.  CLUSTERING 

In  the  previous  section  we  introduced  the  clustering  concept.  In 
this  section  we  shall  provide  more  details  of  the  concept  as  well  as 
the  structure  and  operation  of  various  clustering  models.  Figure  6 
depicts  a  typical  cluster  (in  this  case,  three  trips). 

Clustering  makes  sense  if  the  origins  are  reasonably  close  to¬ 
gether  and  the  destinations  are  also  close  together.  One  way  to 
develop  an  evaluation  of  such  circumstance  is  to  locate  the  centroid 
of  the  origins,  the  centroid  of  the  destinations  and  accumulate  the 
resulting  distances  from  the  original  trips.  In  Figure  6  we  could 
evaluate  the  distances  represented  by  (a+b+c)  +  (d+e+f ) .  If  this 

sum  is  small  then  it  would  make  sense  to  cluster  the  trips. 

By  utilizing  surrogate  distances  we  lose  the  actual  route  dis¬ 
tance  evaluation;  however,  we  gain  the  ability  to  evaluate  large 
numbers  of  cluster  possibilities  conveniently  and  simultaneously. 

In  Figure  6  we  might  employ  Euclidean  distances.  In  this  case  we 
could  compare  the  sum  of  the  row  prices,  p^  +  p^  +  P£,  generated  in 
the  covering  model  for  clustering  to  the  quantity  2(a+b+c)  +  2(d+e+f) 

+  g  to  determine  whether  clustering  is  appropriate.  We  can  think  of 
the  latter  quantity  as  a  surrogate  for  the  vehicle  routing  distance. 

We  can  develop  a  straight  forward  extension  of  the  Location  - 
Allocation  model  discussed  in  Section  7  to  identify  good  clusters 
for  the  dial-a-ride  problem.  In  place  of  the  quantity 
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we  would  employ  a  quantity 


J  l  2[ (x . -a . ) 2  +  (y.-b.)2]  z, . 
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where  x.,  y.,  a.  and  b.  correspond  to  origins  and  x.,  y , ,  and  a. 

.1  J  i  i  j  j  i 

correspond  to  destinations  of  trips. 

It  is  also  possible  to  develop  a  savings  approach  to  clustering 
in  a  similiar  manner  to  that  described  earlier  for  the  delivery  problem. 
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11.  CHAINING 

Upon  termination  of  the  clustering  process,  a  number  of  reasonably 
good  clusters  are  available.  This  set  includes  not  only  the  "best" 
cluster  sets  as  selected  by  the  covering  model  for  clustering,  but 
also  a  number  of  other  clusters  which  might  be  nearly  as  good,  but 
which  were  not  selected  in  the  optimal  solution  to  the  covering  model. 
Figure  7  illustrates  a  set  of  clusters  which  might  result  from  the 
clustering  process.  The  clusters  in  the  figure  represent  only  the 
optimal  solution  to  the  clustering  process.  In  addition,  we  might 
have  other  clusters  available,  e.g.,  a  cluster  containing  only  trips 
1  and  2. 

The  clusters  obtained  represent  good  possibilities  for  segments 
(legs)  of  a  vehicle  route.  The  next  step  in  the  process  is  to  link 
("chain")  these  clusters  into  complete  vehicles  routes.  Figure  8 
illustrates  the  chaining  concept.  Trips  1  and  2  form  one  cluster, 
while  trips  3,  4  and  5  form  another  cluster  (see  Figure  8a).  In 
Figure  8b,  trips  1  and  2  are  replaced  by  a  single  pseudo  (cluster) 
trip,  as  is  also  the  case  for  trips  3,  4  and  5.  These  cluster  trips 
are  then  chained  together. 

The  interpretation  of  chaining  is  that  a  single  vehicle  will 
service  the  first  set  of  trips  (in  Figure  8b  these  would  be  trips  1 
and  2)  and  then  proceed  to  service  the  next  set  of  trips  (i.e.,  trips 
3,  4  and  5)  in  the  chain.  Figure  9  illustrates  the  likely  vehicle 
route  to  service  the  two  clusters. 

We  shall  demonstrate  how  mathematical  models  can  be  developed  to 
aid  the  human  in  developing  good  chains  (vehicle  routes) . 


The  most  fundamental  model  of  chaining  is  based  on  network  flows. 

In  a  network  we  represent  each  cluster  trip  by  two  nodes  (an  origin 
and  a  destination  node)  and  an  arc.  We  then  add  additional  arcs  to  the 
network  which  link  the  destination  node  of  one  cluster  arc  to  the  origin 
node  of  another  cluster  arc  on  the  basis  of  proximity  considerations. 

In  Figure  10,  we  indicate  one  network  flow  model  for  the  clusters  of 
Figure  7. 

To  establish  the  flow  variables  for  the  network  flow  problem  consider 
the  case  where  each  trip  appears  in  exactly  one  cluster.  We  assign  every 
arc  an  upper  capacity  of  1,  the  dashed  arcs  a  lower  capacity  of  0  (zero) 
and  the  solid  cluster  arcs  a  lower  capacity  of  1  indicating  that  these 
trips  must  be  serviced.  Arc  flow  costs  for  the  network  flow  problem 
are  the  associated  vehicle  travel  distances.  If  the  "starting  node" 

and  the  "ending  node"  represents  the  vehicle  storage  depot,  then  we  may 
also  assign  additional  costs  to  the  arcs  originating  at  the  starting 
node  to  reflect  the  fixed  cost  of  using  each  vehicle  (provided  they 
arc  all  t  tie  same )  . 

If  the  network  flow  problem  for  chaining  contains  no  circuits  and 
if  the  same  trip  does  not  appear  in  two  different  clusters  in  the  network, 
then  the  resulting  chains  (vehicles  routes)  will  be  valid  ones.  If  either 
of  these  two  conditions  are  explicitly  included  in  the  model,  then  the 
underlying  network  flow  problem  becomes  considerably  more  complicated. 

We  can  circumvent  some  of  the  problems  associated  with  the  difficult 
problems  of  chaining  by  utilizing  a  partitioning  model  of  chaining.  (Actually 
we  substitute  one  difficulty  for  another.)  In  a  partitioning  model  of  chaining. 
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we  associate  a  column  of  the  partitioning  matrix  for  each  feasible  chain 
(vehicle  route)  and  a  row  for  each  trip.  Column  j  contains  1  in  row 
i  if  trip  i  is  serviced  by  chain  (vehicle  route)  j.  Otherwise,  it 
contains  a  0  (zero).  The  zero-one  variable  associated  with  the 
columns  provide  indications  of  which  chains  were  selected  in  the  optimal 
partitioning  solution.  Table  1  illustrates  a  partitioning  matrix 
for  several  chains  in  Figure  10. 

besides  so  Ivin.',  the  partitioning  model  optimally  (which  we  would  not 
likely  do),  the  obvious  difficulty  with  such  a  model  is  in  generating 
the  candidate  chains  (columns).  To  aide  in  this  column  generation 
process,  we  could  employ  (1)  the  human  (which  would  certainly  be  part 
ot  any  interactive  process) ,  (2)  a  savings  approach  Co  combine  clusters 
into  chains,  (  J)  a  network  flow  model  or  (4)  some  other  method.  The 
first  two  of  these  column  generation  approaches  should  be  reasonably 
intuitive  to  develop.  We  shall  briefly  describe,  the  third  one. 

Suppose  we  have  identified  some  candidate  columns  (chains)  for  the 
partitioning  problem  (e.g.  the  individual  clusters  from  the  clustering 
process).  Then,  the  partitioning  model  will  viel.d  (1)  an  optimal 

» 

solution  and  (2)  a  set  of  row  prices.  As  before,  a  column  j  not  in  the 
solution  appears  favorable  if 

I 
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Consider  how  these  two  terms  might  be  represented 
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Suppose  we  associate  with  each  dashed  arc  the  negat  i  ve  at  wha  t  - 
i'ViT  i-.isis  iii.it  are  incurred  in  traversing  the  arc  (this  is  the  same  as 
before  with  a  sign  change).  Associated  wi Lli  each  clustering  arc 
t  ho  ijuantitv  )  |>  -  c,  where  the  term  }  p  represents  the  sum  of  row 
prices  tor  trips  in  tin-  cluster  and  c  represents  the  cost  of  the  cluster. 
.  i  tli  these  costs  <lef  ined  on  the  network,  we  seel;  a  pain  (or  paths)  in 
the  network,  from  the  starting  to  the  ending  node,  which  has  positive 
total  cost.  Tt  should  he  clear  that  such  a  path  satisfies  the  condition 
for  a  potentially  improving  chain,  and  can  be  added  to  the  partitioning 
proh lorn . 

Again,  with  this  network  flow  model  (a  shortest  path  model  for  a 
sin,  |c  chain)  we  have  the  difficulties  associated  with  circuits  and 
with  the  same  trip  serviced  several  times  by  a  single  chain.  The 
phenomenon  oi  a  trip  being  serviced  several  times  seems  to  occur  in- 
lrec|uently  since  in  seeking  a  least  cost  solution,  Lhe  model  tends  to 
avoid  such  an  instance.  When  it  does  occur,  an  attractive  heuristic 
is  to  simply  delete  one  ol  t  lie  conflicting  clusters  from  the  column 
being  generated.  We  are  currently  utilizing  the  human  interactor  to 
handle  the  case  where  the  flow  problem  contains  circuits.  The  human 
breaks  tin-  circuits  and  patches  the  paths  back  together.  If  we  re- 


strict  the  procedure  to  locating  a  single  improving  path  each  time. 


the  resulting  shortest  path  model  lends  itself  more  intuitively  to 
the  development  of  good  heuristic  procedures  for  solving  both  of  the 
major  difficulties  associated  with  the  flow  mode.  In  particular  it 
is  much  easier  tor  the  human  to  break  the  circuits  in  a  single  path 
than  in  multi paths. 
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An  interactive  dial-a-rido  system  based  on  the  concepts  presented 


here  has  been  implemented  on  a  Chromatics  coLor  graphics  terminal  inter¬ 
faced  with  a  CYBER  74  mainframe  computer.  Although  the  system  is  in 
many  ways  very  rudimentary,  it  dramatically  indicates  the  potential 
for  this  kind  of  interactive  optimization.  We  are  in  the  process  of 
making  extensive  modifications  in  the  software  and  are  testing  a 
variety  of  new  models  and  heuristics  to  aid  in  route  generation. 
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